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1.  Introduction 


Four  basic  methods  to  predict  aerodynamic  forces  and  moments  on  a  projectile  in  atmospheric 
flight  are  commonly  used  in  practice:  empirical  methods,  wind  tunnel  testing,  computational  fluid 
dynamics  (CFD)  simulation,  and  spark  range  testing.  Empirical  methods  have  been  found  very 
useful  in  conceptual  design  of  projectiles  where  rapid  and  inexpensive  estimates  of  aerodynamic 
coefficients  are  needed.  These  techniques  aerodynamically  describe  the  projectile  with  a  set  of 
geometric  properties  (diameter,  number  of  fins,  nose  type,  nose  radius,  etc)  and  catalog  aerody¬ 
namic  coefficients  of  many  different  projectiles  as  a  function  of  these  features.  These  data  are  fit 
to  multivariable  equations  to  create  generic  models  for  aerodynamic  coefficients  as  a  function  of 
these  basic  projectile  geometric  properties.  The  database  of  aerodynamic  coefficients  as  a  function 
of  projectile  features  is  typically  obtained  from  wind  tunnel  or  spark  range  tests.  Examples  of  this 
approach  to  projectile  aerodynamic  coefficient  estimation  include  missile  DATCOM,  PRODAS, 
and  AP981  ( /  through  6).  The  advantage  of  this  technique  is  that  it  is  a  general  method  applicable 
to  any  projectile.  However,  it  is  the  least  accurate  method  of  the  four  methods  mentioned, 
particularly  for  new  configurations  that  fall  outside  the  realm  of  projectiles  used  to  form  the  basic 
aerodynamic  database. 

Wind  tunnel  testing  is  often  used  during  projectile  development  programs  to  converge  on  fine 
details  of  the  aerodynamic  design  of  the  shell  (7,  8).  In  wind  tunnel  testing,  a  specific  projectile 
is  mounted  in  a  wind  tunnel  at  various  angles  of  attack  with  aerodynamic  forces  and  moments 
measured  at  various  Mach  numbers  via  a  sting  balance.  Wind  tunnel  testing  has  the  obvious 
advantage  of  being  based  on  direct  measurement  of  aerodynamic  forces  and  moments  on  the 
projectile.  It  is  also  relatively  easy  to  change  the  wind  tunnel  model  to  allow  detailed  parametric 
effects  to  be  investigated.  The  main  disadvantage  of  wind  tunnel  testing  is  that  it  requires  a  wind 
tunnel  and  is  therefore  modestly  expensive.  Furthermore,  dynamic  derivatives  such  as  pitch  and 
roll  damping  as  well  as  Magnus  force  and  moment  coefficients  are  difficult  to  obtain  in  a  wind 
tunnel  and  require  a  complex  physical  wind  tunnel  model. 

Over  the  past  couple  of  decades,  tremendous  strides  have  been  made  in  the  application  of  CFD  to 
predict  aerodynamic  loads  on  air  vehicles  including  projectiles.  These  methods  are  increasingly 
being  used  throughout  the  weapon  development  cycle  including  early  in  a  program  to  create 
relatively  low  cost  estimates  of  aerodynamic  characteristics  and  later  in  a  program  to  supplement 
and  reduce  expensive  experimental  testing.  In  CFD  simulation,  the  fundamental  fluid  dynamic 
equations  are  numerically  solved  for  a  specific  configuration.  The  most  sophisticated  computer 
codes  are  capable  of  unsteady  time-accurate  computations  with  the  use  of  Navier-Stokes  equations. 
Examples  of  these  tools  include,  for  example,  CFD++,  Fluent,  and  Overflow-D.  CFD  is  computa¬ 
tionally  expensive,  requires  powerful  computers  to  obtain  results  in  a  reasonably  timely  manner, 

*DATCOM  is  not  an  acronym;  PRODAS  =  Projectile  Design  and  Analysis  System;  AP98  =  AeroPrediction98. 
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and  requires  dedicated  engineering  specialists  to  drive  these  tools  (9  through  24).  Spark  range 
aerodynamic  testing  has  long  been  considered  the  gold  standard  for  projectile  aerodynamic 
coefficient  estimation.  It  is  the  most  accurate  method  for  obtaining  aerodynamic  data  on  a  specific 
projectile  configuration.  In  spark  range  aerodynamic  testing,  a  projectile  is  fired  through  an 
enclosed  building.  At  a  discrete  number  of  points  during  the  flight  of  the  projectile  (<  30),  the 
state  of  the  projectile  is  measured  via  spark  shadowgraphs  (25  through  29).  The  projectile  state 
data  are  subsequently  fit  to  a  rigid  six-degree-of-freedom  (6-DOF)  projectile  model  with  the  use  of 
the  aerodynamic  coefficients  as  the  fitting  parameters  (30,  31,  32).  Although  this  technique  is  the 
most  accurate  method  for  obtaining  aerodynamic  data  on  a  specific  projectile  configuration,  it  is 
usually  the  most  expensive  alternative,  requires  a  spark  range  facility,  and  (strictly  speaking)  is 
only  valid  for  the  specific  projectile  configuration  tested. 

Various  researchers  have  used  CFD  to  estimate  aerodynamic  coefficient  estimation  of  projectiles. 
Early  work  focused  on  Euler  solvers  applied  to  steady  flow  problems,  while  more  recent  work  has 
solved  the  Reynolds-averaged  Navier-Stokes  equations  and  large  eddy  simulation  Navier-Stokes 
equations  for  steady  and  unsteady  conditions  (9  through  24).  For  example,  to  predict  pitch 
damping,  Weinacht  prescribed  projectile  motion  to  mimic  a  typical  pitch-damping  wind  tunnel  test 
in  a  CFD  simulation  to  estimate  the  different  components  of  the  pitch-damping  coefficient  of  a  fin- 
stabilized  projectile  (33).  Excellent  agreement  between  computed  and  measured  pitch  damping 
was  attained.  Algorithm  and  computing  advances  have  also  led  to  the  coupling  of  CFD  codes  to 
projectile  rigid  body  dynamics  (RBD)  codes  to  simulate  free  flight  motion  of  a  projectile  in  a  time- 
accurate  manner.  Aerodynamic  forces  and  moments  are  computed  with  the  CFD  solver  while  the 
free  flight  motion  of  the  projectile  is  computed  by  the  integration  of  the  RBD  equations  of  motion. 
Sahu  achieved  excellent  agreement  between  spark  range  measurements  and  a  coupled  CFD-RBD 
approach  for  a  fin-stabilized  projectile  (34).  Projectile  position  and  orientation  at  down-range 
locations  consistent  with  a  spark  range  test  were  extracted  from  the  output  of  the  CFD-RBD 
software  to  compute  aerodynamic  coefficients.  Standard  range  reduction  software  was  used  for 
this  purpose  with  good  agreement  obtained  when  contrasted  against  sample  spark  range  results. 

The  ability  to  accurately  compute  projectile  aerodynamics  in  highly  unsteady  conditions  has  led  to 
the  notion  of  “virtual  wind  tunnels”  and  “virtual  fly-outs”  where  the  simulation  tools  are  used  to 
replicate  a  wind  tunnel  or  spark  range  test. 

Computation  time  for  accurate  coupled  CFD-RBD  simulation  remains  exceedingly  high  and  does 
not  currently  represent  a  practical  method  for  typical  flight  dynamic  analysis  such  as  impact  point 
statistics  (e.g.,  circular  error  probable)  computation  where  thousands  of  “fly-outs”  are  required. 
Furthermore,  this  type  of  analysis  does  not  allow  the  same  level  of  understanding  of  the  inherent 
underlying  dynamics  of  the  system  that  RBD  analysis  with  aerodynamic  coefficients  yields. 
However,  the  coupled  CFD-RBD  approach  does  offer  an  indirect  way  to  rapidly  compute  the 
aerodynamic  coefficients  needed  for  rigid  6-DOF  simulation.  During  a  time-accurate  CFD-RBD 
simulation,  aerodynamic  forces  and  moments  and  the  full  rigid  body  state  vector  of  the  projectile 
are  generated  at  each  time  step  in  the  simulation  (34).  This  means  that  aerodynamic  forces, 


2 


aerodynamic  moments,  position  of  the  mass  center,  body  orientation,  translational  velocity,  and 
angular  velocity  of  the  projectile  are  all  known  at  the  same  time  instant.  With  time-synchronized 
air  load  and  state  vector  information,  the  aero-dynamic  coefficients  can  be  estimated  with  a  simple 
fitting  procedure.  This  report  creates  a  method  to  efficiently  generate  a  complete  aerodynamic 
model  for  a  projectile  in  atmospheric  flight  with  four  short  time  histories  at  each  Mach  number  of 
interest  with  an  industry  standard  time-accurate  CFD-RBD  simulation.  The  technique  is  exercised 
on  sample  CFD-RBD  data  for  a  small  fin-stabilized  projectile.  Parametric  trade  studies 
investigating  the  required  length  of  each  time  snippet  as  well  as  the  required  CFD  accuracy  are 
reported. 


2.  Projectile  CFD-RBD  Simulation 


The  projectile  CFD-RBD  algorithm  employed  here  combines  a  rigid  6-DOF  projectile  flight 
dynamic  model  with  a  three-dimensional  (3-D),  time-accurate  CFD  simulation.  The  RBD  dynamic 
equations  are  integrated  forward  in  time  where  aerodynamic  forces  and  moments  that  drive  the 
motion  of  the  projectile  are  computed  via  the  CFD  algorithm.  The  RBD  projectile  model  allows 
for  three  translational  degrees  of  freedom  and  three  rotational  degrees  of  freedom.  As  shown  in 
figures  1  and  2,  the  I  frame  is  attached  to  the  ground  while  the  B  frame  is  fixed  to  the  projectile 
with  the  IB  axis  pointing  out  the  nose  of  the  projectile  and  the  JB  and  KB  unit  vectors  fonning  a 

right-handed  triad.  The  projectile  state  vector  is  comprised  of  the  inertial  position  components  of 
the  projectile  mass  center  ( x,y,z  ),  the  standard  aerospace  sequence  Euler  angles  the  body 

frame  components  of  the  projectile  mass  center  velocity  (u,v,w),  and  the  body  frame  components 
of  the  projectile  angular  velocity  vector  (p,q,r  ). 


Both  the  translational  and  rotational  dynamic  equations  are  expressed  in  the  projectile  body 
reference  frame.  The  standard  rigid  projectile,  body  frame  equations  of  motion  are  given  by 
equations  1  through  4. 
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Note  that  the  total  applied  force  components  (X,Y,Z)  and  moment  components  ( L,M,N  ) 
contain  contributions  from  weight  and  aerodynamics.  The  aerodynamic  portion  of  the  applied 
loads  in  equations  3  and  4  is  computed  with  the  CFD  simulation  and  passed  to  the  RBD 
simulation. 

On  the  other  hand,  the  CFD  flow  equations  are  integrated  forward  in  time  where  the  motion  of  the 
projectile  that  drives  flow  dynamics  are  computed  with  the  RBD  algorithm.  The  complete  set  of 
3-D  time-dependent  Navier-Stokes  equations  is  solved  in  a  time-accurate  manner  for  simulation  of 
free  flight.  The  commercially  available  code,  CFD++,  is  used  for  the  time-accurate  unsteady  CFD 
simulations  (35,  36).  The  basic  numerical  framework  in  the  code  contains  unified  grid,  unified 
physics,  and  unified  computing  features.  The  3-D,  time-dependent  Reynolds-averaged  Navier- 
Stokes  (RANS)  equations  are  solved  with  the  following  finite  volume  equation. 

—  J  WdV  +  \^(F-G)dA  =  J  HdV  (5) 

v  V 

in  which  W  is  the  vector  of  conservative  variables,  F  and  G  are  the  inviscid  and  viscous  flux 
vectors,  respectively,  H  is  the  vector  of  source  tenns,  V  is  the  cell  volume,  and  A  is  the  surface 
area  of  the  cell  face.  A  second  order  discretization  is  used  for  the  flow  variables  and  the  turbulent 
viscosity  equation.  The  turbulence  closure  is  based  on  topology-parameter-free  formulations. 
Two-equation  higher  order  RANS  turbulence  models  are  used  for  the  computation  of  turbulent 
flows.  These  models  are  ideally  suited  to  unstructured  bookkeeping  and  massively  parallel 
processing  because  of  their  independence  from  constraints  related  to  the  placement  of  boundaries 
and/or  zonal  interfaces. 

A  dual  time-stepping  approach  is  used  to  integrate  the  flow  equations  to  achieve  the  desired  time 
accuracy.  The  first  is  an  “outer”  or  global  (and  physical)  time  step  that  corresponds  to  the  time 
discretization  of  the  physical  time  variation  term.  This  time  step  can  be  chosen  directly  by  the  user 
and  is  typically  set  to  a  value  to  represent  1/100  of  the  period  of  oscillation  expected  or  forced  in 
the  transient  flow.  It  is  also  applied  to  every  cell  and  is  not  spatially  varying.  An  artificial  or 
“inner”  or  “local”  time  variation  term  is  added  to  the  basic  physical  equations.  This  time  step  and 
corresponding  “inner  iteration”  strategy  is  chosen  to  help  satisfy  the  physical  transient  equations 
to  the  desired  degree.  For  the  inner  iterations,  the  time  step  is  allowed  to  vary  spatially.  Also, 
relaxation  with  multi-grid  (algebraic)  acceleration  is  employed  to  reduce  the  residues  of  the 
physical  transient  equations.  It  is  found  that  an  order  of  magnitude  reduction  in  the  residues  is 
usually  sufficient  to  produce  a  good  transient  iteration. 

The  projectile  in  the  coupled  CFD-RBD  simulation,  along  with  its  grid,  moves  and  rotates  as  the 
projectile  flies  down  range.  Grid  velocity  is  assigned  to  each  mesh  point.  For  a  spinning  and 
yawing  projectile,  the  grid  speeds  are  assigned  as  if  the  grid  were  attached  to  the  projectile  and 
spinning  and  yawing  with  it. 
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In  order  to  properly  initialize  the  CFD  simulation,  two  modes  of  operation  for  the  CFD  code  are 
used,  namely,  an  uncoupled  and  a  coupled  mode.  The  uncoupled  mode  is  used  to  initialize  the 
CFD  flow  solution  while  the  coupled  mode  represents  the  final  time  accurate  coupled  CFD-RBD 
solution.  In  the  uncoupled  mode,  the  RBD  are  specified.  The  uncoupled  mode  begins  with  a 
computation  performed  in  “steady  state  mode”  with  the  grid  velocities  prescribed  to  account  for 
the  proper  initial  position  ( x0 ,  y0 ,  z0 ),  orientation  ( (f>t] ,  00 ,  t//0 ),  and  translational  velocity  ( u0 ,  v0 ,  w0  ) 

components  of  the  complete  set  of  initial  conditions  to  be  prescribed.  After  the  steady  state 
solution  is  converged,  the  initial  spin  rate  ( )  is  included  and  a  new  quasi-steady  state  solution 

is  obtained  with  time-accurate  CFD.  A  sufficient  number  of  time  steps  are  performed  so  that  the 
angular  orientation  for  the  spin  axis  corresponds  to  the  prescribed  initial  conditions.  This  quasi¬ 
steady  state  flow  solution  is  the  starting  point  for  the  time-accurate  coupled  solution.  For  the 
coupled  solution,  the  mesh  is  translated  back  to  the  desired  initial  position  (x0,y0,z0)  and  the 

remaining  angular  velocity  initial  conditions  ( q0,r0 )  are  then  added.  In  the  coupled  mode,  the 

aerodynamic  forces  and  moments  are  passed  to  the  RBD  simulation  which  propagates  the  rigid 
state  of  the  projectile  forward  in  time. 


3.  Flight  Dynamic  Projectile  Aerodynamic  Model 


The  applied  loads  in  equations  3  and  4  contain  contributions  from  projectile  weight  and  body 
aerodynamic  forces  and  moments: 
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The  terms  containing  CWA  constitute  the  Magnus  air  load  acting  at  the  Magnus  center  of  pressure 
while  the  terms  containing  CX0,CX2,CNA  define  the  steady  load  acting  at  the  center  of  pressure. 

The  externally  applied  moment  about  the  projectile  mass  center  is  composed  of  an  unsteady 
aerodynamic  moment  along  with  terms  because  the  center  of  pressure  and  center  of  Magnus  are 
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not  located  at  the  mass  center.  The  tenns  involving  CMA  account  for  the  center  of  pressure  being 
located  off  the  mass  center  while  the  terms  involving  CNPA  account  for  the  center  of  Magnus  being 

located  off  the  mass  center.  The  aerodynamic  coefficients  are  all  a  function  of  local  Mach  number 
which  are  typically  handled  through  a  table  look-up  scheme  in  projectile  flight  simulation  codes. 


4.  Projectile  Aerodynamic  Coefficient  Estimation  (PACE) 


The  time-accurate  coupled  CFD-RBD  simulation  provides  a  full  flow  solution,  including  the 
aerodynamic  portion  of  the  total  applied  force  and  moment  (X,Y,Z,L,M,N)  along  with  the  full 
state  of  the  rigid  projectile  ( x,  y,  z,  </>,  6 ,  y/,  u,v,  w,  p ,  q ,  r  )  at  every  time  step  in  the  solution  for  each 
time  snippet.  Given  a  set  of  n  short  time  histories  (snippets)  that  each  contain  m  time  points,  we 
obtain  a  total  of  h  =  m  *  n  time  history  data  points  for  use  in  estimating  the  aerodynamic  coef¬ 
ficients:  Cx0 ,  Cx 2 ,  CNA ,  CYPA  ,  CLDD ,  CLP ,  CMA ,  CMQ ,  CNPA .  Note  that  for  fin-stabilized  proj  ectile 

configurations,  the  Magnus  force  and  moment  are  usually  sufficiently  small  so  that  CWA  and  CNPA 
are  set  to  zero  and  removed  from  the  fitting  procedure  to  be  described  next. 

Equations  6  and  7  represent  the  applied  air  loads  on  the  projectile  expressed  in  the  projectile  body 
frame.  Computation  of  the  aerodynamic  coefficients  is  aided  by  transfonnation  of  these  equations 
to  the  instantaneous  aerodynamic  angle  of  attack  reference  frame  that  rotates  the  projectile  body 
frame  about  the  I B  axis  by  the  angle  ;/  =  tan1  (  h7  v  )  . 
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(9) 


Each  time  history  data  point  provides  a  total  of  six  equations  given  by  the  components  of  equations 
8  and  9.  The  first  component  of  equation  8  is  gathered  together  for  all  time  history  data  points  to 
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form  equation  10.  Likewise,  the  second  and  third  components  of  equation  8  generate  equations  1 1 
and  12,  respectively,  while  the  first  component  of  equation  9  constructs  equation  13.  Finally,  the 
second  and  third  components  of  equation  9  are  gathered  together  to  fonn  equation  14.  Subscripts 
on  the  projectile  state  vector  and  aerodynamics  force  and  moment  components  represent  the  time 
history  data  point. 
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Equations  10  through  14  represent  a  set  of  five  uncoupled  problems  to  solve  for  the  different  aero¬ 
dynamic  coefficients.  To  estimate  the  aerodynamic  coefficients  near  a  particular  Mach  number,  a 
set  of  n  time  accurate  coupled  CFD-RBD  simulations  are  created  over  a  relatively  short  time 
period.  Since  an  individual  time  snippet  is  over  a  short  time  period  where  the  projectile  state 
variables  does  not  change  appreciably,  it  is  critical  that  initial  conditions  for  the  different  time 
snippet  be  selected  in  an  infonned  way  so  that  the  rank  of  each  of  the  fitting  matrices  is  maximal. 
Properties  of  the  fitting  matrices,  such  as  the  rank  or  condition  number,  can  be  used  as  an  indicator 
of  the  suitability  of  the  CFD-RBD  simulation  data  to  estimate  the  aerodynamic  coefficients  at  the 
target  Mach  number.  Equation  10  is  employed  to  estimate  the  zero  yaw  drag  coefficient  ( Cxo)  and 

the  yaw  drag  coefficient  ( CX2 ).  To  minimize  the  condition  number  of  this  fitting  matrix,  both  low 

and  high  aerodynamic  angle  of  attack  time  snippets  are  required.  Equation  1 1  is  used  to  compute 
the  normal  force  coefficient  ( CNA  )  and  it  requires  time  history  data  with  a  non-zero  aerodynamic 

angle  of  attack.  Equation  12  is  used  to  compute  the  Magnus  force  coefficient  ( CWA )  and  it  requires 

time  history  data  with  both  low  and  high  roll  rate  and  aerodynamic  angle  of  attack.  Equation  13  is 
employed  to  estimate  the  fin  cant  roll  coefficient  ( CLDD  )  along  with  the  roll-damping  coefficient 

( CLP  ).  To  minimize  the  condition  number  of  this  fitting  matrix,  both  low  and  high  roll  rate  time 

snippets  are  required.  Equation  14  is  employed  to  estimate  the  pitching  moment  coefficient 
( CMA  ),  the  pitch-damping  coefficient  ( CMQ  ),  and  the  Magnus  moment  coefficient  ( CNPA  ).  For 

successful  estimation  of  these  coefficients,  time  history  data  with  both  low  and  high  roll  rate  and 
aerodynamic  angle  of  attack  as  well  as  low  and  high  aerodynamic  angle  of  attack  are  required.  To 
meet  all  the  requirements  for  successful  estimation  of  all  five  sets  of  aerodynamic  coefficients,  four 
time  snippets  are  used,  all  with  different  initial  conditions.  Table  1  lists  the  four  cases  with  launch 
conditions.  Notice  that  the  set  of  time  snippets  contains  a  diverse  set  of  initial  conditions:  zero 
aerodynamic  angle  of  attack  and  angular  rates;  high  angle  of  attack  and  zero  angular  rates;  low 
angle  of  attack,  high  roll  rate  with  other  angular  rates  zero;  zero  angle  of  attack,  high  pitch  rate  with 
other  angular  rates  zero. 
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Table  1.  Time  snippet  initial  conditions. 


State 

Case  1 

Case  2 

Case  3 

Case  4 

x(m) 

0 

0 

0 

0 

y  (m) 

0 

0 

0 

0 

z(m) 

0 

0 

0 

0 

Phi  (deg) 

0 

0 

0 

0 

Theta  (deg) 

0 

0 

0 

0 

Psi  (deg) 

0 

0 

0 

0 

V(m/s) 

1032 

1032 

1032 

1032 

v  (m/s) 

0 

0 

0 

0 

w  (m/s 

0 

-352.5 

-90 

0 

p  (rad/s) 

0 

0 

377 

0 

q  (rad/s) 

0 

0 

0 

-10 

r  (rad/s) 

0 

0 

0 

0 

Alpha  (deg) 

0 

20 

5 

0 

For  flight  dynamic  simulation,  aerodynamic  coefficients  are  required  at  a  set  of  Mach  numbers 
that  covers  the  intended  spectrum  of  flight  conditions  for  the  round.  If  aerodynamic  coefficients 
are  estimated  at  k  different  Mach  numbers,  then  a  total  of  l  =  k*n  time  snippets  must  be 
generated  to  construct  the  entire  aerodynamic  database  for  flight  simulation  purposes. 


5.  Results 


In  order  to  exercise  the  method  developed,  a  generic  finned  projectile  is  considered.  A  sketch  of 
the  projectile  is  shown  in  figure  3.  The  projectile  has  the  following  geometric  and  mass  properties: 
length  =  0.1259  m,  reference  diameter  =  0.013194  m,  mass  =  0.0484  kg,  mass  center  location  from 
base  =  0.0686  m,  roll  inertia  =  0.74E-06  kg  m2,  pitch  inertia  =  0.484E-04  kg  m2. 


Figures  4  through  9  present  projectile  state  trajectories  for  each  of  the  four  time  snippets.  Each 
time  snippet  is  0.023  second  and  contains  50  points,  leading  to  an  average  output  time  step  of 
0.0004.  The  initial  conditions  for  each  of  the  time  snippets  are  shown  in  table  1.  These  four 
snippets  create  time  history  data  at  low  and  high  angle  of  attack,  roll  rate,  and  pitch  rate  needed  for 
accurate  aerodynamic  coefficient  estimation.  Notice  that  Cases  2  and  3  have  notably  more  drag 
because  of  the  high  angle  of  attack  launch  conditions.  Case  3  is  launched  with  relatively  high  roll 
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rate  compared  to  all  other  cases.  Case  4  generates  roll  rate  toward  the  end  of  the  time  snippet 
because  of  high  angle  of  attack  roll-pitch  coupling.  Significant  oscillations  in  Euler  pitch  angle  are 
created  in  Case  2  with  some  cross-coupling  response  exhibited  in  Euler  yaw  angle.  Figures  10 
through  1 5  plot  aerodynamic  forces  and  moments  in  the  local  angle  of  attack  reference  frame 
defined  for  Cases  1,3,  and  4  since  these  cases  are  the  primary  ones  used  to  estimate  the  coef¬ 
ficients.  For  all  cases,  the  axial  force  oscillates  from  -20  N  to  -25  N.  There  exists  a  slight  bias 
between  the  CFD-RBD  and  estimated  data  of  about  0.5  N  for  low  angle  of  attack  time  snippets. 

For  moderately  high  angles  of  attack  (Case  3),  the  estimated  data  also  oscillate  with  a  much  higher 
amplitude  than  the  CFD-RBD  data,  indicating  that  CX2  is  estimated  larger  than  the  CFD-RBD 
suggests.  The  normal  force  time  snippets  agree  well  between  the  CFD-RBD  and  estimated  data 
for  all  time  snippets.  For  the  sample  finned  projectile,  side  force  and  out-of-plane  moment  are 
generally  small  (<  0.5  N,  0.05  Newton  meters  [Nm])  because  of  a  negligibly  small  Magnus  force 
and  moment.  The  CFD-RBD  and  estimated  data  agree  reasonably  well  but  certainly  do  not  over- 
lie  one  another.  The  only  time  snippet  that  creates  notable  rolling  moment  is  Case  3  which  is 
launched  with  an  initial  roll  rate  of  377  rad/sec.  Notice  that  the  estimated  data  smoothly  go 
through  the  CFD-RBD  data  which  oscillate  in  a  slightly  erratic  manner.  The  in-plane  moment 
(Mz)  agrees  reasonably  well  for  both  the  CFD-RBD  and  estimated  data.  The  results  shown  in 
figures  4  through  15  are  typical  for  all  Mach  numbers.  The  overall  observation  from  the  data  is 
that  the  estimated  aerodynamic  model  fits  the  CFD-RBD  data  well,  with  the  notable  exception  of 
axial  force  where  a  bias  between  the  two  is  exhibited. 


Figure  4.  Velocity  for  the  time  snippets. 
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Figure  8.  Euler  pitch  angle  for  the  time  snippets. 
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Figure  9.  Euler  yaw  angle  for  the  time  snippets. 


-15 


-20 


-25 


-30 


Cases  1  a  4,  CFD/RBD 
/ 

Case  1  a  4,  Estimated 

1 

. „M;;..>iuvg;i|||||'lH 

. <MMIM||||||^|II|U|llll 

\  Case 

,  CFD/RBD 

Case  3,  Estimated 

_ 1 _ 

0  0.005  0.01  0.015  0.02  0.025 

Time  (sec) 


Figure  10.  Estimated  (dashed)  and  CFD-RBD  (solid)  body  axis  axial  force  (Fx) 
versus  time. 
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Figure  12.  Estimated  (dashed)  and  CFD-RBD  (solid)  side  force  (Fz)  versus  time. 
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Figure  13.  Estimated  (dashed)  and  CFD-RBD  (solid)  body  axis  rolling  moment  (Mx) 
versus  time. 
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Figure  14.  Estimated  (dashed)  and  CFD-RBD  (solid)  pitching  moment  (My)  versus  time. 
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The  sample  projectile  investigated  in  this  report  has  been  fired  in  a  spark  range  at  Mach  3.0  with 
aerodynamic  coefficients  computed  via  conventional  aerodynamic  range  reduction.  Table  2 
presents  a  comparison  of  aerodynamic  coefficients  obtained  from  spark  range  testing  and  sub¬ 
sequent  coefficients  obtained  with  the  method  described  here.  Notice  that  most  aerodynamic 
coefficients  such  as  Cxo,  Cna,  and  Cma  are  in  reasonably  good  agreement  with  the  test  data.  Axial 
force  yaw  drag  and  roll  damping  are  both  different  by  -20%,  while  pitch  damping  is  different  by 
-40%.  Differences  in  the  spark  range  and  estimated  coefficients  can  be  attributed  to  several 
factors,  such  as  inaccuracies  in  the  CFD-RBD  solution,  the  manner  in  which  the  spark  range 
data  were  reduced,  and  inaccurate  estimation  because  of  insufficiently  rich  data. 


Table  2.  Comparison  of  estimated  aerodynamic  coefficients  and  estimated  coefficients  at  Mach  3.0. 


Spark  Range 
Data  -  Spark 
Range  Reduction 

CFD-RBD  - 
PACE 

Percent  Difference 
Between  Coefficients 

Zero  Yaw  Axial  Force  Coefficient,  Cxo 

0.221 

0.238 

7.1  to  7.7 

Yaw  Axial  Force  Coefficient,  Cx2 

5.0 

5.9 

15.0  to  18.0 

Nonnal  Force  Coefficient  Derivative, 

Cna 

5.83 

5.64 

3.2  to  3.3 

Pitching  Moment  Coefficient  Derivative, 

Cma 

-12.6 

-13.82 

8.8  to  9.7 

Pitch-Damping  Moment  Coefficient,  Cmo 

-196 

-134 

31.6  to  46.3 

Roll-Damping  Moment  Coefficient,  CLP 

-2.71 

-3.37 

19.6  to  24.4 

CFD-RBD  data  were  generated  at  six  different  Mach  numbers  ranging  from  1.5  to  4.0.  The 
estimation  algorithm  discussed  was  used  to  compute  a  complete  set  of  aerodynamic  coefficients 
across  its  Mach  range.  These  results  are  provided  in  table  3  with  plots  of  the  individual  aerody- 
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namic  coefficients  given  in  figures  16  through  21.  With  the  exception  of  Cxi,  the  steady  aero¬ 
dynamic  coefficients  are  smooth  and  follow  typical  trends  for  variation  in  Mach  number.  The  yaw 
drag  coefficient,  Cx2,  however,  is  somewhat  erratic  with  a  low  value  of  0.21  at  Mach  1.5  followed 
by  a  steady  rise  until  Mach  3.5.  Pitch  damping  decreases  with  Mach  number,  as  would  be 
expected  for  a  fin-stabilized  projectile  beyond  Mach  1.0.  However,  roll  damping  steadily  increases 
until  Mach  4.0  when  in  drops  off  notably. 


Figure  16.  Zero  yaw  axial  force  coefficient  versus  Mach  number. 


Figure  17.  Yaw  axial  force  coefficient  versus  Mach  number. 
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In  order  to  investigate  convergence  of  the  estimation  procedure,  the  number  of  points  used  per  time 
snippet  was  varied  from  a  single  point  per  time  snippet  to  50  points  per  time  snippet.  Of  all  the 
coefficients,  the  pitch-damping  coefficient  generally  required  the  most  number  of  points  to 
converge.  Figure  22  plots  estimated  pitch-damping  coefficient  versus  the  number  of  points  per 
time  snippet.  For  this  aerodynamic  coefficient,  convergence  is  reached  after  about  50  points  per 
time  snippet.  The  points  extracted  from  each  time  snippet  are  from  the  beginning  of  the  time 
history  and  are  not  evenly  extracted  over  the  entire  time  history. 


Figure  22.  Pitch-damping  moment  coefficient  versus  number  of  points  per  time  snippet. 


6.  Conclusions 


With  a  time-accurate  CFD  simulation  that  is  tightly  coupled  to  an  RBD  simulation,  a  method  to 
efficiently  generate  a  complete  aerodynamic  description  for  projectile  flight  dynamic  modeling  is 
described.  A  set  of  n  short  time  snippets  of  simulated  projectile  motion  at  m  different  Mach 
numbers  is  computed  and  employed  as  baseline  data.  The  combined  CFD-RBD  analysis  computes 
time-synchronized  air  loads  and  projectile  state  vector  infonnation,  leading  to  a  straightforward 
fitting  procedure  to  obtain  the  aerodynamic  coefficients.  The  estimation  procedure  decouples  into 
five  sub-problems  that  are  each  solved  via  linear  least  squares.  The  method  has  been  applied  to  a 
sample  supersonic  finned  projectile.  Overall,  the  results  are  encouraging.  A  comparison  of  spark 
range  obtained  aerodynamic  coefficients  with  the  estimation  method  presented  here  at  Mach  3 
exhibits  good  agreement  within  10%  for  Cxo,  CNA,  and  CMa;  agreement  within  20%  for  CX2  and 
Clp;  and  poor  agreement  within  40%  for  Cmq.  Convergence  of  the  aerodynamic  coefficients  is  a 


21 


strong  function  of  the  number  of  points  in  each  time  snippet  with  pitch  damping  generally  requiring 
the  most  number  of  points  for  convergence.  This  technique  reported  here  provides  a  promising 
new  means  for  the  CFD  analyst  to  predict  aerodynamic  coefficients  for  flight  dynamic  simulation 
purposes.  It  can  easily  be  extended  to  flight  dynamic  modeling  of  different  control  effectors, 
provided  that  accurate  CFD-RBD  time  simulation  is  possible  and  an  aerodynamic  coefficient 
expansion  is  defined  which  includes  the  effect  of  the  control  mechanism. 
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